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TECHNICAL NOTE D- 


LONGITUDINAL OSCILLATION OF 
LAUNCH VEHICLES 


INTRODUCTION 


Any rocket (missile or launch vehicle) represents an elastic structure. During flight, any 
disturbance caused by wind gusts or irregularities of the thrust may cause longitudinal oscillations of 
the structure. Iii general oscillations of this kind are of a transient nature and are damped out after 
a certain time. The vibrating system is said to be stable. 

It is well known, however, that under certain circumstances the system proves to be un- 
stable; that is, vibration energy can be tapped from the thrust by the activity of the system itself and 
sustain the structural oscillations. The amplitudes increase until some nonlinear effects limit any 
further increase. In other words, self-excited oscillations, so called “POGO” oscillations named for 
the pogo-stick-like motion of the vehicle, may result. 

The first time, POGO oscillations occurred during test flights of the Titan-II missile [1] . 

Later, the first and second stages of the Saturn V Launch vehicle were plagued by POGO oscillations [2] . 

The supply of energy is caused by an interaction of the structure and the propellant system. 
Oscillation of the structure, especially of the pumps, causes pressure oscillation of the liquid propellants. 
Consequently, an oscillating pump discharge and hence an oscillating thrust is generated that in return 
can excite the structure. In this way, the feedback loop is closed (Fig. 1). An inherent mechanism of 
the system may cause feedback of the thrust oscillation to sustain the structural oscillation at least for a 
certain time interval (Fig. 2). 

To examine the stability of the vehicle, 
a proper model must be established. Then the 
characteristic equation must be investigated to 
determine whether there are roots having positive 
real parts that imply instability. Commonly this 
will be done by application of an appropriate 
stability criterion. In this way, the stability 
conditions’ in terms of the system parameter can 
be obtained. 


Since the POGO problem developed, 
numerous stability analyses have been performed. 
A review of the efforts made for the Saturn 
V/ Apollo vehicle is given in Reference 3. In 
References 3 and 4 a description of one of the 
latest models and its stability analysis is out- 
lined. According to these reports the POGO 
Figure 1. Feedback system with constant time lag t. model is based on the idea that a model must 

be sufficiently complex to serve as a prediction 

tool. It should cover all experimental and flight data. In other words, the basic idea is that the more 
accurately the model copies the actual system, the better it will display all the important features, 
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Figure 2. S-IC POGO oscillations; SA-502; Apollo 6. 


especially the one of stability. Accordingly, the model is a coupled system of structural modes, propellant 
line modes, and inboard and outboard engine gain and phase transfer functions. The structural modes 
were obtained from a 280-degree-of-freedom structural model. The propellant line modes were derived 
also from many degree-of- freedom systems composed essentially of lumped liquid masses connected by 
springs. The stability analysis leads to a set of equations that forms a closed-loop system. The roots of 
the characteristic equations were obtained by the root locus technique. 

Although opinions about the validity of such complex models may differ, one peculiarity 
is evident. Even if one assumes the model and analysis to be correct, one cannot consider both of them 
to be completely satisfying. The reason is that no understanding of the mechanism and the actual 
physical cause of self-excitation can be obtained in this manner. Any decision concerning the stability of 
the system is left to the computer. 

In this note, the POGO problem shall be examined from another point of view. Two special 
features of the system prove to be crucial for creation of instability. 

1. There is an inherent time lag of the thrust oscillations behind the structural oscillations 
(Fig. 1). The time lag arises from the (constant) time interval required to convert the pump inlet pressure 
oscillations into thrust oscillations. 

2. The decrease of mass caused by the exhausting of gas. 1 


1. During the burn phase a vehicle loses approximately 70 percent of its mass. 
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Systems with time lag or systems possessing retarded action, as they are sometimes referred 
to, are known and have been examined previously [5-8] . The point is that systems with time lag are 
capable of self-excited oscillations or instability. In general, retarded action is present in all closed-loop 
control systems [5,8] . However, an example of a mechanical system was given by Minorsky in 1942 
[5] . The equations of motion of problems possessing retarded action are amenable to so-called 
difference-differential equations [5,7] . This notation stems from the fact that these equations contain 
at least one term, called a retarded term, which depends not on the time t itself but on the difference 
t-r, where r represents the time lag. In its linear form, having one degree-of-freedom only, the problem 
is already of transcendental character. It always leads to an infinite spectrum of frequencies with which 
the system can oscillate. In practice however, only a few of them (most frequently only one) are actually 
responsible for self-excitation [7]. 

To illustrate the importance of the time lag and difference-differential equation concepts for 
systems possessing retarded actions, Minorsky may be quoted. In essence he pointed out [5] : 

“In fact in all closed-loop systems there always exists a time lag in 
the cyclic process of operation. A mathematical description of such a process 
in terms of a difference-differential equation simplifies the commonly estab- 
lished procedure of ‘couplings’ between the different parts of the system which 
generally leads to an ordinary differential equation of a relatively high order; 
this usually involves considerable difficulties if one attempts to discuss the 
characteristic equation algebraically. Furthermore the nature of coupling 
very frequently remains unknown. On the other hand, in establishing a 
difference-differential equation of a retarded process one must know only 
the over-all time lag r, without the necessity of knowing the intermediate 
links in the system.” 


The conclusions to be drawn for the problem at hand are obvious. The combined work of 
the single parts of the propellant system must not be considered. Only the system’s total effect on the 
vibrating structure is of significance. This effect is determined by two quantities, however: the time lag 
and the magnitude of the feedback. As far as the structure is concerned, a number of structural modes 
may serve for its representation. 

The model derived in this note follows this line. However, the model is intended to serve 
first for understanding the basic mechanism of self-excitation; thus, further simplifications have been 
made. Accordingly, nonessential system details that may efface the actual cause of instability have been 
omitted. In doing so, the model has been reduced to one structural mode. This is in accordance with 
experience because all flight measurements show one POGO frequency only. Another simplification 
concerns the oscillating thrust component. Although the entire propellant system occurs in its creation, 
only the oscillation of the pumps is considered. The importance of the latter is shown in Reference 9. 

By linear approximation the respective amplitudes are assumed to be proportional. The thrust oscillation, 
however, is delayed by the time interval r. 

Accordingly, the simple model proposed is a single mass-spring model possessing retarded 
action. It has been obtained by a one-term Galerkin approach to the equation of motion. The stability 
is examined with the aid of Nyquist’s criterion [6,8] . Special emphasis has been placed on showing the 
effect of the mass decrease. The possible importance of the latter was quoted in Reference 10. Further- 
more, it is shown that a significant conclusion can be drawn from this simple model. 
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At this point, an important fact must be mentioned. The actual vibrating system is a non- 
linear one whereby its nonlinearity represents an essential feature. The model derived thust obviously be 
a linear one which covers a small vibration of the system only. As is well known, the linear theory is in- 
sufficient to describe the actual motion of the system. Especially, limit cycles cannot be determined in 
this way. By using linear analysis, the conditions only can be obtained under which the system is stable 
or unstable. 


EQUATION OF MOTION 


During powered flight the vehicle mass center is moving along the flight path. If an oscillating 
thrust component AT is acting on the vehicle, it causes vibrational motion of the vehicle mass center 
according to Newton’s first law. Consequently, inertia forces of the vehicle masses are created that equi- 
librate the oscillating thrust component AT. This equilibrium system of forces, however, excites free 
longitudinal vibration of the elastic vehicle structure that superimposes the rigid body motion. If the 
vehicle structure is reduced to a lumped mass-spring system, the equation of motion of the vibrating masses 
relative to the mass center is 


Mx(t) + Dx(t) + Kx(t) = ATe n - x Q (t) MX q 


a) 


where, according to Newton’s first law, 


AT = M 0 x 0 and M 0 


n 



( 2 ) 


Substitution of equation (2) into equation (1) results in 


Mx + Dx + Kx = AT 




(3) 


Equation (3) shall be solved by a Galerkin approach [1 1]. Restriction to one term yields the 
desired one-degree-of-freedom model mentioned previously. 

The eigensolutions of the above problem, 


coq — 0 ; Xq 


1 

1 

1 


(rigid body motion) 
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and 


coj #0 ; Xj 


x,0)‘ 


X 2 


0) 


0) 


; j = l,2,...n-l 


(4) 


are the solutions of the matrix eigenvalue problem 


(K - co 2 M) x = 0 


(5) 


The eigenvectors satisfy the orthogonality conditions: 


Xj'MXj- J m k x k (i) x k G) = SjjMj ; i,j = 0,1,2, ...n-1 , 


( 6 ) 


where M 0 has been defined by equation (2). From equations (5) and (6) it follows that 


Xj' KX- = ojj 2 5“Mj = 5jj Kj , i,j = 0,l,2,...n-l 


(7a) 


hence, 




co:” = ; j = 0,l,2,...n-l 

1 M; 


(7b) 


To perform the Galerkin approach, it is assumed that 


n 1 

x(t) = aj(t) Xj 


( 8 ) 


where aj (j= s l,...n-l) has the dimension of a length. Substitution of equation (8) into equation (3) 
and left multiplication by the eigenvector of the eigenvectors (4) yields 


n-1 . 


M k a k + “2 aj X k ' DXj + K fc a k = XTxJ® ; k = l,2,...n-l ; (9) 

thereby, equations (6) and (7a), are considered. 
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In Reference 12, it is explained that no considerable error results if one assumes 


X^. DXj — > kj - 1,2, ...n-1 


Thus, it follows from equation (7b) and (9) 

x n (k ) 

a k (t) + 2e k a k (t) + ,o> k 2 a k (t) = — — AT ; k = l,2,...n-l (10) 


Now, AT must be expressed. As mentioned previously, it is assumed that the oscillating 
thrust is caused primarily by the pump oscillations. In addition it is assumed that the magnitudes of 
both oscillations are linear-dependent. Because AT represents the retarded quantity that lags the time 
interval r behind the structural oscillations, one may write 


AT = - sx^t - t) 


(ID 


where x fl is the coordinate of the mass at station n to which the pumps belong and the factor of 
proportionality s has the dimension of a mass. 

The- acceleration of the pumps can be composed of that of the mass center and relative 
acceleration. The latter is given by equation (8). Thus, 


n-1 


= x r 


+ 2 




a) 


( 12 ) 


However, the magnitude of the oscillating thrust AT is not large enough to accelerate the vehicle mass 
center considerably. In other words, x c is a small quantity compared with the relative acceleration and 
can therefore be neglected. ' 

Then, from equations (1 1) and (12), it follows that 


AT = - s 


n-1 

2 

1=1 


aj(t - r) XjG) 


(13) 


Substitution of equation (13) into equations (10) yields the following system of equations: 

V*) + 2e k + w k 2 a k^) = -JT x n (k) 2 a/t-r) x n C) ; k = 1,2 ... n-1 

k j=l 


Because a one-term approach is considered, the system reduces to 


6 



a'i(t) + 2e 1 a 1 (t) + co^a^t) = - ^ x^ 1 ) a\(t-r) 


(14) 


According to equations (6) and (7a) 


M, 


x or 

n . 



K, 



(15) 


where m t * , k x * represent generalized mass and spring constant of the first mode normalized to the unit 

at station n. One realizes that the smaller x is relative to the other coordinates of the first mode, 
the larger is the generalizes mass m? : 

If the subscript 1 is omitted, one obtains from equations (14) and (15) 

la(t) + 2fcoa(t) + co 2 a(t) = -r?a(t-r) (16a) 


thereby are 


? = 



£l 

CO 


V 55 




(16b) 


(16c) 


where 



(16d) 


In the following, the stability of the simple system represented by equations (16) will be examined. 

APPLICATION OF NYQUISTS STABILITY CRITERION 


By the assumption 
a(t) = ae^ 

it follows from equation (16a) that 


(17a) 
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0 


(17b) 


X 2 + 2ftoA + to 2 + X 2 17 = 


This is a transcendental equation having an infinite number of roots that, in general, are complex: 


\j - jtf® ± i7® ; j = 1,2... 


(17c) 


If interpreted in the complex A-plane, one may state: The vibration is stable [6, 8] if and only if no 
lies in the right half-plane. In other words, the vibration is stable if 


< 0 ; j = 1,2... 

Because the spring-mass system which follows from the system (16a) by omitting the right 

side 


a + 2fcoa + co 2 a = 0 


can be assumed to be a vibrating system with positive (not overcritical) damping, its characteristic 
equation 


X 2 + 2?c oX + co 2 = 0 


(18) 


does not have roots in the right half-plane. Consequently, one may divide equation (17b) by the left side 
of the equation (18) without gaining any poles in the right half-plane. Thus, instead of equation (17b), 
the following equation may be examined: 


1 +p 


X 2 e- Xr 

X 2 + 2f&)X + co 2 


0 


or 


1 + f(X) = 0 


where 


f(X) = V 


X 2 e~ Xr 

X 2 + 2?coX + bj 2 ' 


(19) 
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To study the stability of the system at hand, the Nyquist analysis must be applied [8]. The 
Nyquist criterion is: 

The system is stable if and only if f(X) does not encircle the 
point (-1,0) of the f(X)-plane if X travels along the so-called Nyquist path. 

The Nyquist path is a closed contour in the X-plane that completely encloses the entire right 
half of the X-plane. In Reference 6 it has been shown that steady-state solutions of equation (16a) exist 
under very special conditions only. Thus, the general case, 


M® * 0 ; j = 1,2,- (20) 


shall first be studied. In doing so, one may consider the Nyquist path to be going along the imaginary 
axis [6, 8]. To enclose the right half-plane, a large semicircular path is drawn in the right half-plane having 
the radius R that is interpreted as being infinite in the limit (Fig. 3). 

Because 


lim X 2 e = 0 
R *> 


and, hence, 


lim f(X) = 0 
R 00 

Figure 3. Nyquist path in the 
complex X-plane. 

the traveling along the semicircle does not contribute 
(if R -» 00 ) and therefore can be omitted [6]. 

Thus, if one substitutes 


X = icon (21a) 


the Nyquist path is represented by the following: 


-00 < a < co 


From equations (19) and (21a) one realizes that 

2 -iojra 

f(ic off) = 7 ? (21b) 

a 2 - 1 - 2i fa 
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Now f(icoa) shall be expressed by its radius vector r(o) and its argument. First, the 
denominator: 


a 2 - l - 2$o = V ( a 2 - l) 2 + 4f 2 a 2 e“^ 


where 


iT*t> = 


cos 0 - i sin 0 


sin 0 - 


V (a 2 - l) 2 + 4f 2 a 2 


and 


cotan 0 = 


J_ a 2 - 1 
2f a 


Thus, equation (21b) can be written 


f(icoa) = r(a) “ ojra ) 


where 


r(u) = i? 


V(o 2 - l) 2 + 4f 2 a 2 


and 0 is given by equations (22) and (23). 
Because 


( 22 ) 


(23) 


(24a) 


(24b) 


f(-icoa) = f(icoa) 


whereby the conjugate complex f is symmetrical to f, the real axis being the symmetry axis, one may 
restrict the travel path to positive a. Thus, the Nyquist criterion can be expressed as follows: 
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The system is stable if and only if f(A) does not encircle the point (-1,0) of the f(A)-plane if 
a travels along ' 


0 < o < (25) 

where X and o are linked by equation (21a). 

Now, one has to determine under what conditions f(X) encircles '(-1,0) if a travels along 
the interval (25). 

Because 


1 = e ±i (2j - 1> 


j = 1,2 ... 


it follows from equations (24) that encirclement will occur if and only if a positive o exists so that 

<t> - uto = ±(2j - l)7r ; j = 1,2 ... (26) 

where 0 is determined by equations (22) and (23), and 


r(a) > 1 

By combining equations (22), (23), (25), and (26), it follows: 

, 1 a 2 - 1 

cotan cora = — 

2 $ a 

(2j - 1 )n < cora < 2j7r ; j = 1 , 2 ,... 

Equation (28) is a transcendental equation having infinite many roots. Let 


(27) 


(28) 

(29) 


Oj j = 1,2,3,.. 


(30) 


be those positive roots which belong to the intervals (29). Then, considering equation (24b) and the con- 
dition (27), one may state: Encirclement will occur if for at least one root a k of the roots (30) 


V 2 > 




(31a) 
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No encirclement can occur if for all roots (30) 


V < 1 


2r 


CT j 


j = 1,2,3- 


(31b) 


Figure 4 is a graph of 


gw 


■ ("*)■ •(?) 


(32) 


The function g(a) represents the right side of the conditions (31). It shows two stationary points at 


m 


■S- 


1 -2f 2 


and 


The respective values of g(a) are 


g(tf m > = d-f 2 ) 


g(oo) = 1 

The first represents a minimum, and the second represents a maximum. The horizontal line through 
o = 1 is the asymtote of g(°°). 
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The abscissa a j , o ^ are the positive solutions of 


g(a) = n 


(33a) 


One obtains 



1 - 2£ 2 
1 -V 2 


1 ± 



l - n 

(1 - 2 f 2 ) 2 


(33b) 


For any Oy. satisfying 


a T < a, < 
I k 


'II 


(34) 


condition (31a) is valid as can be seen from Figure 3. Thus, encirclement can occur if and only if the 
interval (34) contains at least one of the roots (30). 

Let n Q be the number of encirclements of (- 1 , 0); then, from Figure 4 or equation 
(33b), one concludes 


n = 0 
o 


if r? < 2 ? sf\ - T 2 


(35a) 


n = 00 if r] > 1 

o 


(35b) 


In the latter case, namely, the length of the interval (34) becomes infinite. 


The conclusions which can be drawn from Figure 4 and equation (33b) are gathered in 

Table 1. 


Now, the special case of the steady-state oscillation shall be examined. It occurs if at least 
one of the real parts (20) equals zero. Then, according to equation (17c), one may assume 


X = ± iy ; 7 > 0 


(36) 


Now, from equations (16), (17), and (36), it follows that 


a + ( 2 £cu + r?7sin7r)a + (co 2 - T77 2 cos7T)a = 0 . 
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For this equation to represent the steady-state oscillation, the coefficients must reduce to 
2 f co 4* 177 sin 77 = 0 
co 2 - 177 2 cos 7 r = 7 2 
or if notation (21a) is used ( o > 0 ) 

7] o sin co t o = -2f (37a) 


77 o cos co 7 o 



Equations (37) are equivalent to the following equations: 


cotan co 7 a 




(37b) 
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( 38 ) 


(2 j - 1) 7 r < COT <J < 2 j 7T 

■ (S^)’ + (-£)’ 


These equations, in turn, are identical to equations (28), (29), (32), and (33a). Thus, the 
special conditions under which steady-state oscillations may occur are characterized by coincidence of 
one of the roots (30) with either o j or In that case, 


7i = wor ! or y 2 = 


are the frequencies of steady-state oscillations. 

Now, combining this and previous statements, especially equations (35), Nyquist’s criterion 
for the case at hand is: 

1. The system is stable if 


V < 2? Vl ~ 2 

2. If 

T? > 2$Vi -r 2 


the system is stable if none of the roots (30) lie within the interval 


a 


I 




a < a 


II 


(39) 


which is determined by equations (33). If at least one of the roots (30) belongs to this interval, the 
system is unstable. 

3. If rj > 1, the system is unstable. 

For some practical applications, another formulation of Nyquist’s criterion is appropriate. It follows 
from equation (22), (26) conditions of encirclement (31), and equation (38) that the system is 
stable or unstable accordingly as 


n < — -2L - — » 
a sm co r o 
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for all of the roots (30) or 


-M — 

asm cora 


for at least one of these roots. 

Figure 4 is a graphic interpretation of the stability problem. For determination of the roots 
(30), the curves (28) are drawn in. Here, cotan cora is represented by the first and the second branch 
of the intervals (29), j = 1, 2. Figure 4a shows stability because o x lies to the right of a^ . Conse- 
quently all roots (30) lie outside the interval (30). However, Figure 4b is an unstable case because o x 
lies between a^ and a^ . 



(a) STABLE CASE 



Figure 5. Stability graphs. 
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STABILITY OWING TO THE DECREASE IN MASS 


One realizes that the stability of the system depends, apart from the damping ratio f , on 
two dimensionless parameters, namely 

cor and o o r 
s 

While the first one determines the roots (30), the ratio 



determines the boundaries o ■ and of the critical interval (39). 

The quantities s and r represent specific quantities of the vibrating system and can 
therefore be considered to be constant. Now, the decrease in mass may cause a change over from the 
state of stability (Fig. 4a) to one of instability (Fig. 4b) and vice versa. The reason is that a decrease 
in mass affects both of the parameters r? and co r. The decrease in mass causes, in general, an increase 
of the natural period c o [ 1 1 ] 2 and a decrease of the generalized mass m* as it can be seen by a 
closer examination . 

At the beginning of the powered flight, the liquid propellant masses are large; consequently, 
77 is a small quantity. The increase of 17 may first cause o ^ and to become real and later on 

may increase the length - a j of the interval (39) as shown in Figure 3 and Table 1. As result of 

the increase of co , however, the roots (30) (Fig. 4) travel in the negative direction of the a-axis. 
Consequently one of the roots, in general the first one, approaches the interval (39) and finally may 
enter it. However, this implies instability. Further decrease of mass may cause this root to leave the 
interval. If at that time the next one is still outside the interval, the system regains stability; otherwise, 
the instability persists. In this way, the state of stability may change over to one of instability and vice 
versa. 

To demonstrate this behavior a numerical case shall be considered. Assumed is a single 
mass-spring system (m* = m). The damping ratio of the system shall be 


? = 0.06 


The state of decreasing mass shall be characterized by an interval of the natural frequency co 
4.6 < cor < 5.9 (40) 


2. In exceptional cases, co remains constant as the mass decreases. 
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The Nyquist criterion is applied to some successive (discrete) values of cot from this interval. Five 
cases determined by different values of cor are investigated. All values of cot and co r are 
gathered in Table 2. s s 


TABLE 2. STABILITY STATE OF A ONE-DOF SYSTEM DEPENDENT ON NATURAL 
FREQUENCY AND FEEDBACK MAGNITUDE. 


COT 

COT 

S 

1 

2 

3 

4 

5 

14.40 

14.33 

14.25 

14.00 

13.40 

4.60 

S 

s 

S 

s 

S 

5.00 

S 

s 

S 

u 

u 

5.17 

S 

s 

u 

u 

u 

5.30 

S 

u 

U 

u 

u 

5.39 

S 

s 

u 

u 

u 

5.60 

S 

s 

s 

u 

u 

5.90 

S 

s 

s 

s 

u 


Note: S = stable; U = unstable. 


According to the stability criterion, if 
cot 

co t > / = 

s V 1 " f 


the system is stable. Thus, for the above f and the upper limit of the interval (40) (cor = 5.9), it 
follows that for 


co r > 17.34 
s 


the system is stable for the entire interval (40). Because the co s r values of Table 2 are smaller than 

the above value, special consideration must be given to the position of the roots (30) relative to the 
respective intervals (39) according to the stability criterion. The result of this investigation is shown by 
Table 2. 

Cases 2 through 5 of Table 2 show transition from stability to instability and, with 
exception of the last one, return to stability. 


1'8 




The Nyquist criterion is not an algebraic criterion. It is based on fundamentals of the theory 
of the complex variables. For that reason, it lacks some perspicuity. It therefore seems to be con- 
venient to perform also a direct integration of equation (16a). This has been done for the cases 3 and 
4 of Table 2. 

The integration is performed by the method of Runge-Kutta (Table 3 and 4) and also by 
analog computer (Figs. 6 and 7). It is assumed that 


r = 1. 

The initial conditions are 


a(o) = .2 ; a(o) = 0. 


As far as the Runge-Kutta analysis is concerned, the amplitude maxima only are printed (Table 3 and 4). 


TABLE 3. CASE 3 (co $ = 14.25) MAXIMA OF a(t). 


w 
















5.0 

0.200 

0.151 

0.152 

0.153 

0.151 

0.151 

0.150 

0.150 

0.149 

0.149 

0.148 

0.148 

0.147 

0.147 

0.145 

53 

0.200 

0.148 

0.153 

0.153 

0.154 

0.155 

0.155 

0.156 

0.157 

0.157 

0.157 

0.158 

0.158 

0.159 

0.159 

5.6 

0.200 j 

0.144 

0.150 

0.149 

0.148 

0.147 

0.146 

0.146 

0.145 

0.145 

0.143 

0.143 

0.142 

0.142 

0.141 


Owing to the transient condition at the beginning, the values of a(t) behave nonuniform, but, 
later, the tendency to stability or instability, respectively, can easily be observed. 

TABLE 4. CASE 4 (co § = 14.00) MAXIMA OF a(t). 


CO 
















4.6 

0.200 

0.152 

0.150 

0.145 

0.141 

0.135 

0.131 

0.128 

0.125 

0.121 

0.117 

0.113 

0.110 

0.107 

0.103 

5.3 

0.200 

0.149 

0.156 

0.157 

0.160 

0.162 

0.165 

0.167 

0.168 

0.171 

0.174 

0.176 

0.180 

0.182 

0.183 

5.9 

0.200 

0.140 

0.146 

0.141 

0.140 

0.137 ; 

0.134 

0.132 

0.129 

0.128 

0.125 

0.124 

0.121 

0.118 

0.117 


Once again, one realizes the transition from stability to instability and the return to stability, 
caused by the decrease of mass. However, these results are obtained by linear analysis which, as men- 
tioned previously, is insufficient to describe the actual motion. As in the case of Nyquist’s criterion, 
the conditions only are displayed under which the system is stable or unstable. 
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SUMMARY AND CONCLUSION 


Because of its elastic structure, a launch vehicle is capable of longitudinal oscillation. In 
this regard, three basic features prove to be essential: 

1. The elastic vehicle structure and the propellant system form a feedback system. 

2. A time delay exists in the feedback loop. 

3. The loss of mass because of the depletion of liquid propellant. 

The intent of this report is to show that the presence of these features impart the tendency to self- 
excitation. 

It is well known that feedback systems with time lag are sensitive to self-excitation or 
instability. However, the loss of mass changes the systems mass and mass distribution continuously. 
Consequently, the entire system changes continuously and also may change its stability behavior. At 
one time point, the system may be stable, and at another time point, it may be unstable. 

From this point of view, the changing of the stability state of a vehicle during powered 
flight must be considered a certainly peculiar but absolutely understandable feature of the vehicle. It 
may be difficult or impossible to design a vehicle which proves to be stable during the entire powered 
flight time. Above all, it seems to be very difficult to predict the stability behavior of a vehicle. The 
only remaining alternative is to keep the effect of self-excitation unsignificant. Preventive measures 
exist. 


First, the quantity rj should remain small; consequently, the effect of feedback and also 
the interval of instability (39) remain small. This can be achieved by large m*. According to the 
remark following equation (15), the coordinate x n must be kept as small as possible relative to the other 

coordinates of the first mode. This implies to support the engines as stiff as possible. The connection 
between soft engine support and self-excitation was demonstrated in the past. First, the connection was 
demonstrated during test flights of the Titan-II missile, whereby POGO oscillations first occurred. 

The spring support of the engine-pump masses of this missile is striking [1] . Later, POGO oscillation 
of the S-II stage during some of the Apollo flights were observed. Especially, severe oscillations occurred 
during the flight of Apollo 13 [2] . Flight data obtained at that time show extremely large relative 
displacements of the center engine. The soft support of the center engine by a pin-ended cross beam 
led to this situation. 

The use of gas-filled accumulators [1,3] in some of the suction lines represents another 
preventive measure. An accumulator is a device which is able to absorb pressure fluctuations in the 
suction system before they could propagate through the pumps and produce significant thrust 
oscillations. In the case of the Titan-II missile and also of the S-IC stage of the Saturn V launch 
vehicle, these “POGO fixes” proved to be successful [1, 3] . 
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